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Abstract 



We study the phase transition of a system of self-gravitating neutrinos in the pres- 
ence of a large radiation density background in the framework of the Thomas-Fermi 
model. We show that, by cooling a non-degenerate gas of massive neutrinos below 
some critical temperature, a condensed phase emerges, consisting of quasi-degenerate 
supermassive neutrino stars. These compact dark objects could play an important 
role in structure formation in this universe, as they might in fact provide the seeds 
for galactic nuclei and quasi-stellar objects. 
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A gas of massive fermions which interact only gravitationally has interesting ther- 
mal properties which may have important consequences for the early universe. The 
canonical and grand canonical ensembles for such a system have been shown to have 
a nontrivial thermodynamical limit |2| . Under certain conditions these systems will 
undergo a phase transition that is accompanied by gravitational collapse ||. This 
phase transition occurs uniquely for the attractive gravitational interaction of neutral 
fermions. There is indeed no such phase transition in the case of charged fermions [Q. 
Of course, gravitational condensation will also take place if the fermions have an addi- 
tional weak interaction, as neutrinos, neutralinos, axions and other weakly interacting 
massive particles generally do. Moreover, it is important to note that the phase tran- 
sition will occur, irrespective of the magnitude of the initial density fluctuations in 
the fermion gas, and, as we shall shortly demonstrate, irrespective of the amount of 
background radiation. To be specific, we henceforth assume that this neutral fermion 
is the heaviest neutrino v T , although this is not essential for most of the subsequent 
discussion. 

The ground state of a gravitationally condensed neutrino cloud, with mass below 
the Chandrasekhar limit, is a cold neutrino star ||, [7], ||, in which the degeneracy 
pressure balances the gravitational attraction of the neutrinos. Degenerate stars of 
neutrinos in the mass range between m u = 10 and 25 keV are particularly interest- 
ing H , as they could explain, without resorting to the black hole hypothesis, at least 
some of the features that are observed around supermassive compact dark objects, 
which are reported to exist at the centers of a number of galaxies || [TO], [II], |T~2|, IB, 13 1 



including our own [jT5|, [16|] and quasi-stellar objects (QSO) |T7|, [L8|, |19|, j2"0"| . Indeed, 



there is little difference between a supermassive black hole and a neutrino star of the 
same mass near the Chandrasekhar limit, a few Schwarzschild radii away from the 
object. 

The existence of a quasi-stable neutrino in this mass range is neither ruled out by 
particle and nuclear physics experiments nor by direct astrophysical observations ||. 
In the early universe, however, it would lead to an early neutrino matter dominated 
phase some time after nucleosynthesis and prior to recombination. In such a uni- 
verse, the microwave background temperature would be reached much too early to 
accommodate the oldest stars in globular clusters, cosmochronology and the Hubble 
expansion age, if the Standard Model of Cosmology is correct. However, the early 
universe might have evolved quite differently in the presence of such a heavy neu- 
trino. In particular, it is conceivable that primordial neutrino stars have been formed 
in local condensation processes during a gravitational phase transition that must have 
occurred some time between nucleosynthesis and recombination. Aside from reheat- 
ing the gaseous phase of heavy neutrinos, the latent heat produced by the condensed 
phase might have contributed partly to reheating the radiation as well. Moreover, 
the bulk part of the heavy neutrinos (and antineutrinos) will have annihilated effi- 
ciently into light neutrinos via the Z° in the interior of these supermassive neutrino 
stars ||, ||, 0, §. Since both these processes will increase the age of the universe, or 
the time when the universe reaches today's microwave background temperature j5|, 
it does not seem excluded that a quasi-stable massive neutrino in the mass range 
between 10 and 25 keV is compatible with the cosmological observations |5|, ||. 
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The purpose of this paper is to study the formation of such a neutrino star during a 
gravitational phase transition in an expanding universe at the time when the energy 
densities of neutrino matter and radiation are of comparable magnitude. 

We assume here that the equilibrium distribution of the v T gas is spherically sym- 
metric and the energy density p 7 of the radiation background homogeneous. At this 
stage p 7 consists of photons and the two remaining relativistic neutrino species 
and v ei is given by 

with 



7 / 4 \ 4/3 
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The gravitational potential V(r) satisfies the Poisson equation 

AV = 4:7iGm u (m u n L , + p 7 ), (3) 

where the number density of r neutrinos (including antineutrinos) of mass m v can 
be expressed in terms of the Fermi-Dirac distribution at a finite temperature T as 

njr) = ^^Tf'nl^my (4) 

with 

im - rgu. (5) 

J o 1 + ^ 

g v denotes the combined spin degeneracy factors of neutrinos and antineutrinos (i.e. 
g v is 2 or 4 for Majorana or Dirac neutrinos respectively), and p is the chemical 
potential. It is convenient to introduce the normalized reduced potential 

' G^-n (6) 



m u GM & 

Mq being the solar mass, and dimensionless variable x = t/Rq with the scale factor 

\ ^/^ / \ 8/3 
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EQ-@ then takes the simple form 

where we have introduced the normalized inverse temperature defined as 

P = T /T- T = m u GM Q /R . (9) 
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At zero temperature we recover from eq.(^) the well-known Lane-Emden differential 
equation [6,8] 

d 2 v v 3/2 

% - "71' < 10 > 

The solution of the differential equation (|8|) requires boundary conditions. We assume 
here that the neutrino gas is enclosed in a spherical cavity of radius R corresponding 
to X\ = R/Rq. We further require the total neutrino mass to be M v , and the total 
radiation mass within the cavity to be M 7 , and we allow for the possibility of a 
pointlike mass Mc at the origin, which could be e.g. a compact seed of other exotic 
matter. v(x) and v(x) is then related to its derivative at x = X\ by 

v(xi) = —\v{xi) — 1 , 11 

xi\ M Q J 

which in turn is related to the chemical potential by p = T v'(xi). v(x) at x = is 
related to the point mass at the origin by M c /M & = v(0). 

Similar to the case of the Lane-Emden equation, it is easy to show that eq.(|8|) has 
a scaling property: if v(x) is a solution of eq.@ at a temperature T and a cavity 
radius R, then v(x) = A 3 v(Ax) with (A > 0) is also a solution at the temperatures 
T = v4 4 T, T 7 = v4 4 T 7 and the cavity radius R = R/A. 

It is important to note that only those solutions that minimize the free energy are 
physical. The free energy functional is defined as |J 

F[n] = fi[n]N v -W[n] 

^ r d 3 rd 3 p , /A ( p 2 V[n] M[n]V 

- r Ww M1+exp -2Sj- t + t P- (12) 



where 



and 



(27T) 

3 / m u n(r')+p 1 



V[n] = - Gmv Jd 3 r' m ^^ P \ (13) 
W [n\ = —Gml J d z rd 3 r' U ^^ . (14) 



The chemical potential in eq. (12) varies with density so that the number of neutrinos 
N v = M u /m u is kept fixed. 

All the relevant thermodynamical quantities such as number density, pressure, free 
energy, energy and entropy can be expressed in terms of v/x 



M 3 

nJx) = 



171 u Rq 8tt 



P'^h (^) , (15) 
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F = \»{N V + N,) + *-M 2 1 ± 

+ \T,RlJ d 3 x(n u - n,) v{x) ~ v{0) -RlJ d 3 xP„(x), (17) 

- \t q RI f d 3 x (n l/ + n 1 )^- + (n v -n 1 )^- + R 3 f d 3 xe kin (x), (18) 

S = ^(E-F). (19) 



In eqs. (|T7|) and (|1|) we have introduced the effective radiation number iV 7 = M 1 /m u , 
and the effective radiation number density n 7 = p 1 /m u . 

We now turn to the numerical study of a system of self-gravitating massive neutrinos 
with arbitrarily chosen total mass M = 1OM varying the cavity radius R and the 
radiation mass M 7 . Owing to the scaling properties, the system may be rescaled to 
any physically interesting mass. For definiteness, the v T mass is chosen as m v = 17.2 
keV which is about the central value of the mass region between 10 and 25 keV f§ 
that is interesting for our scenario. 

In fig. 1 we present our results for a gas of neutrinos in a cavity of radius R = 100i?o, 
and with no background radiation i. e. with M 7 = 0. We find the three distinct 
solutions in the temperature interval T = (0.049 — 0.31 l)To of which only two are 
physical, namely those for which the free-energy assumes a minimum. The density 
distributions, corresponding to such two solutions are shown in the first plot in fig. 
1. We refer to the solution that exists above the mentioned interval at arbitrary 
high temperature as "gas", while the solution which exists at low temperatures and 
eventually becomes a degenerate Fermi gas at T = 0, we refer to as "condensate". 

In fig. 1 we also plot various extensive thermodynamical quantities (per neutrino) 
as functions of neutrino temperature. The phase transition takes place at the point 
T t where the free energy of the gas and condensate become equal. The transition 
temperature T t = 0.19442T is indicated by the dotted line in the free energy plot. 
The top dashed curve in the same plot corresponds to the unphysical solution. At 
T = T t the energy and the entropy have a discontinuity. 

Next we study a system of massive neutrino gas in a radiation background. In the 
course of universe expansion, the heavy neutrino becomes nonrelativistic at a time t^R 
corresponding to the temperature T^r = m v . At that time the radiation dominates 
the matter by about a factor of 15. A neutrino cloud with the total mass e.g. 1O 9 M , 
which is by about factor of 10 below the Chandrasekhar limit for m u = 17.2 keV and 
g u = 4, would have a radius 

R NR = 2[G(M„ + M 7 )4 fl ] 1 / 3 , (20) 

yielding Rnr = 0.265 light days. As we shall shortly see this is way below the critical 
value and therefore there exists only one solution for any fixed temperature. As the 
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universe expands and cools down the relative amount of radiation mass decreases as 



R 



M 1 = 15M U ^. (21) 
R 

Concerning the R dependence of the temperature we assume that during the radiation 
domination the temperature T of the neutrino gas is fixed by the radiation heat bath 
so that it continues to decrease asl/R until the neutrino matter starts to dominate. 
At that point the system being in gaseous phase will have the entropy per particle 
S/N v = 7.60. From then on, the T will not be coupled to the radiation heat bath 
anymore and will decrease with R according to the adiabatic expansion of the neutrino 
matter itself. 

Shortly after the matter starts to dominate, we reach the region of instability and 
the first order phase transition takes place. Fig 2. shows how the entropy behaves 
around the critical region. We find as we decrease the radius that the first order 
phase transition becomes weaker and eventually disappears at R = R c = 5.34 light 
days corresponding to the critical temperature T c = 0.0043m,, and the radiation mass 
relative to the neutrino mass M 7 = 0.67M U . This is the critical point of a second 
order phase transition. This behavior is typical for a mean-field type of models [f2~l[ 



The initial entropy of 7.6 per neutrino drops to 0.75 so that the latent heat per particle 
yields AE = TAS = 0.01677V Thus, the condensate formation is accompanied by a 
release of considerable amount of energy which will reheat the radiation environment. 
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Figure Captions: 

1: Density distribution normalized to unity for condensate-like (solid line) and gas- 
like (dashed line) solutions at T — T t . Free energy, energy and entropy per 
particle as a function of temperature. Temperature, energy and free energy, are 
in units of T . 

2: Entropy per particle as a function of temperature for various radii (in light days) 
of the neutrino cloud with M v = 10 9 M Q . 
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Fig. 1 




Fig. 2 



